Spectral interferometric polarised coherent anti-Stokes Raman spectroscopy 
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We have developed an interferometric implementation of coherent anti-Stokes Raman scattering 
(CARS) which enables broadband coherent Raman spectroscopy free from non-resonant background 
(NRB), with a signal strength proportional to concentration. Spectra encode mode symmetry 
information into the amplitude response which can be directly compared to polarised spontaneous 
Raman spectra. The method requires only passive polarisation optics and is suitable for a wide 
range of laser linewidths and pulse durations. 



Spontaneous Raman scattering provides a powerful op- 
tical route to obtain chemically specific information. Its 
application to microscopy enables imaging using molec- 
ular vibrations as a contrast mechanism and mapping of 
cell and tissue constituents based on their chemical signa- 
ture Coherent Raman scattering (CRS) is the non- 
linear multiphoton equivalent and allows much faster ac- 
quisition, [3j, |4| with intrinsic optical sectioning in imag- 
ing However, while CRS imaging has been very suc- 
cessfully employed for imaging individual Raman bands, 
it struggles to harness the powerful chemical specificity 
provided by spectral detection in the Raman fingerprint 
region. In particular, quantitative broadband coherent 
Raman spectroscopy of biological samples, analogous to 
spontaneous Raman spectroscopy, has proven difficult to 
achieve primarily due to the coherent backgrounds in- 
herent to CRS. In this letter we report on a new method 
for quantitative broadband CRS spectral imaging, which 
uses passive polarisation optics combined with spectrally- 
resolved balanced homodyne detection. The technique 
has relaxed requirements on spectral phase and instru- 
ment stability and provides full access to the Raman 
fingerprint region, while retaining the advantages of en- 
hanced signal and optical sectioning inherent to CRS. 

The two most widely used CRS techniques are stim- 
ulated Raman scattering (SRS) and CARS (Fig Q^a)) 
At intensities suitable for biological samples, SRS 
requires heterodyne methods [7|, |8| to detect the sig- 
nal against the coherent background of the excitation 
fields; such methods require wavelength scanning to build 
up spectral information. CARS is a four wave mixing 
(FWM) process with a signal field generated at the anti- 
Stokes frequency ui a s = w 



pr 



— u)s which is spectrally 
separated from the excitation fields (u> p , ujs, i^ pr are de- 
noted the pump, Stokes and probe fields, respectively). 
CARS can therefore simultaneously generate an entire 
vibrational spectrum by using a spectrally broad Stokes 
beam Q, making CARS suitable for rapid detection of 
vibrational spectra, particularly in microscopy. 

The coherent background in CARS arises from FWM 
processes that are independent of vibrational transitions 
(FigQJa)). The anti-Stokes intensity is determined by the 



third order susceptibility, which (away from electronic 
resonances) is given by 



2 XiY H R-e{x fl } + |x R | 



(1) 



where xr an d Xnr are components resonant and non- 
resonant with vibrational modes, respectively. Since 
Xr is proportional to the number of resonant modes 
in a medium, the concentration dependence of CARS is 
quadratic and non-linearly mixed with the non-resonant 
background (NRB). The interferometric (second) term, 
however, is linear in xr ancl is amplified by xnr- Re{xfi} 
is dispersive and antisymmetric about the vibrational 
line-centre; however, Imjx^} is directly related to the 
spontaneous Raman spectrum 

Many methods have been developed to remove the 
effect of the NRB from CARS spectra Non- 
interferometric techniques, such as polarisation-based 
[n]-[l3j and time-resolved methods [3) 15 1 , recover a 



(typically small) proportion of the energy in the third 
term of ([I]). However, the interferometric term is the 
largest source of signal for weak resonances and, below 
the damage threshold of biological samples, CRS has 
been shown to be faster than spontaneous Raman spec- 
troscopy only if this term is detected 0,13. 

Imjxj?} may be recovered through interference be- 
tween the anti-Stokes field and a local oscillator (LO), 
which may be generated either externallv|l7H20|. or 
by using the NRB as an internal LO[21 -23]. Broad- 
band external approaches acquire an entire spectrum 
simultaneously via spectral interferometry, requiring 
high stability[24f and broadband transform limited laser 
sources [2^ (precluding the use of fields generated in pho- 
tonic crystal fibres (PCF), which are typically much 
broader). Using the NRB as the LO ensures that the 
phase relationship between the LO and anti-Stokes fields 
is fixed, which relaxes stability and spectral coherence re- 
quirements; however, as the LO is in phase with Re{x}, 
recovering Im{%} is not trivial. Previous experimen- 
tal implementations also need broad transform limited 
pulses [H, HfiHill , or require phase scanning 21. 28|. 

Computational internal LO approaches have also been 
developed, which use the Xnr ^ e {x R } term present in 
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therefore be used to generate broad, NRB-free, vibra- 
tional spectra. Furthermore, the NRB is removed with- 
out requiring an independent measurement of its spectral 
variation. The technique is similar to Dual Quadrature 



(a) 



Figure 1: (colour online) (a) Energy level diagram. CARS in- 
volves the path uj p — > ljs u-pr —> u a s, SRS is oj p — > u>s — ► 
(where is a vibrational resonace) . An example non-resonant 
FWM pathway is shown as uj p — s> cj' pr — s> oj's — > uj a s- (b) Po- 
larisations for SIP-CARS with one elliptical (or circular) po- 
larised field E e and one linear field Ei . Detection is performed 
along n+ and n_ . E e is shown as the input to a A / 4 plate (ori- 
entated at angle <f> to the fast axis) to generate the desired 
elliptical polarisation. For a circular Stokes/linear pump: 
E e — Es, Ei = E p — E pr , <j> = w /4. For elliptical pump/linear 
Stokes: E e = E p = E pr , Ei = Es, <f> = 7r / sor3 ' r / 8 - ( c ) Ex- 
perimental setup. Both pump/probe (P/Pr: 1050 nm) and 
Stokes (St: 1070-1600 nm) beams were derived from a com- 
mercial supercontinuum source with two outputs (repetition 
rate 20MHz, ~5ps pulses at the sample). Beams were com- 
bined on a dichroic mirror (DM) and focussed into the sample 
by a 1.3NA oil immersion DIC objective lens (Obi), where the 
power of each beam was 30 mW. Samples were prepared be- 
tween two coverslips, separated by <60 fim. The signal was 
collected by a 0.75NA air objective (Ob2), separated from the 
excitation fields by a shortpass filter (SP), and then dispersed 
via grating spectrometer onto a deep depletion CCD. The 
pump/probe ellipticity was controlled by a combination of lin- 
ear polariser (LP) and zero-order half- and quarter-waveplates 
( A /2 and A / 4 > respectively). A Wollaston prism (WP) was 
placed after Ob2 to separate the two detection polarisations 
on the CCD. The beams were again circularly polarised be- 
fore the spectrometer by an achromatic quarter-wave plate 
(not shown) to account for the the different grating reflectiv- 
ity for each polarisation. 



raw CARS spectra to calculate Im{x n } [29. |30|. These 
produce good approximations to Im{x_R} if the spectra 
are of sufficient width, however both methods intrinsi- 
cally produce a spectrally varying error signal, which in- 
creases at resonances (up to 10%) (3lJ. Both techniques 
also require prior or post estimation of the spectral vari- 
ation of the NRB (effectively the spectral variation of 
the Stokes field) from a reference material, which is often 
compromised by the presence of a residual resonant re- 
sponse. Changes of the background spectrum during an 
acquisition leads to errors that can mask the weaker res- 
onances in the fingerprint region of biological samples 

We describe here a new internal LO technique, Spectral 
Interferometric Polarised CARS (SIP-CARS), that is sig- 
nificantly simpler to implement than previous experimen- 
tal approaches. It does not have stringent requirements 
on the lasers used, and is suitable for narrowband, mul- 
tiline and broadband systems. Transform limited pulses 
are not required and PCF generated supercontinua can 



2J, |28|, [321, except 



Spectral Interferometry (DQSI)[2£ 
that the fields have different frequencies and hence only 
interfere via the non-linear response. As a result, the 
two quadratures measured in SIP-CARS (i.e. the real 
and imaginary components of the non-linear response) 
contain different linear combinations of the tensor el- 
ements of x, due to polarisation mixing. The third- 
order response is solved exactly, without assumptions 
on the relative strength of resonant and non-resonant 
components (22, 33 1. 

To illustrate the method, consider CARS with a (right- 
hand) circularly polarised Stokes, and pump and probe 
linearly polarised along the rr-axis, as in Figure [TJb). For 
convenience, we express the Stokes electric field in terms 
of the equivalent linear polarisation, Es , used to gener- 
ate it via a A / 4 plate; so, the field along x- and y-axes 
is Es x = Es/V% and E Sy = iE s /V^ (along any two 
orthogonal directions the Stokes electric field will differ 
in phase by Y 2 )- The susceptibility can be separated 
into diagonal (e.g. Xuu) and off-diagonal (e.g. 
ements; the off-diagonal terms mediate the coupling of 
orthogonally polarised excitations into the detected po- 
larisation. Considering the polarisation, P + , induced in 
the medium along the n + axis, at to the pump and 
probe polarisation, we have 



P+ = lx^E pr E p (E+)* + ~ Xl22X E pr E p (E+y 
+ ^Xii22 E pr E p (Eg)* + -X 1212 EprE p (Eg)* 

where E p and E pr are the electric fields of the pump 
and probe, and Eg and Eg are the components of the 
Stokes along n + and n , respectively (Fig. Dlh)). As 
Eg = iEg , the last two terms lag the first two in phase by 
1T /2. Imaginary components of the last two terms there- 
fore interfere with real components of the first two, with 
the strength of the interference determined by the relative 
strength of the diagonal and off-diagonal tensor elements. 
Similarly, for P_ the Eg and Eg terms are swapped, and 
the last two terms lead the first two by 7T /2. Imaginary 
components are therefore added to the real components 
along n + , and subtracted from them along n_; spectral 
interferometric detection is performed by taking the dif- 
ference between spectra measured at these polarisations , 
leaving the purely imaginary components [22j. 

The induced polarisations are more succinctly ex- 
pressed within a basis including the pump and probe 
polarisations, 



Px — ^XiinEprEpEg 
Py — — 7/2X2 ii^Epr EpEg 



(2) 
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Along the detection axes n + and n_ the induced po- 
larisations are P + = A=(P x +P y ) and P_ = -^(P x — P v ), 
and the anti-Stokes signals are S+ cx P+P+, S_ oc 
p_p* xhe sum and difference of the anti-Stokes in- 
tensities are then £S = S+ + S_ oc P X P X + PyPy an d 
AS = S+ - S_ cx P X P* + (P x Py)*, respectively. The 
difference signal is therefore given by 

AS cx Im {XuuCJWs ( 3 ) 

where li — EiE* are the beam intensities. Separat- 
ing the susceptibilities into resonant and nonresonant 
components x ijk i — X<^ + X^ jkl > assuming an isotropic 

medium (X2112 X12217 Xim X1212 d~ X1221 X1122)' 
and noting that the non-resonant terms possess Klein- 
man symmetry 0| = Xnr, X™ = X™ = X** = 
Xnj?/3) this becomes 

AS cx x„ R Im { x f xii - 3xf 221 } Ws 
cx(l-3p)x„Im{xf m }VV S 

where p = X R 221 Ix^ is the CARS depolarisation ratio 
of the resonance [35]. The difference spectrum is there- 
fore linear in the imaginary component of \, is amplified 
by Xjvhi ancl contains no real, dispersive terms or non- 
resonant contributions. Due to the linear response, well 
established linear multivariate analyses such as principal 
component analysis can be applied. Mode symmetry in- 
formation is mixed into the amplitude response through 
the depolarisation ratio, p (0 < p < 3 / 4 )j and the spec- 
trum can be compared directly to polarised spontaneous 
Raman spectra. 

Because the interference is effectively between 
phase shifted copies of the same fields there are no extra 
requirements on the coherence of the excitation pulses. 
Moreover, if spectra are measured simultaneously, in- 
coherent backgrounds (such as two-photon fluorescence) 
and any variation of the real components of the CARS 
signal are common-mode in S+ and S— and are automat- 
ically subtracted out (though they will still contribute 
shot noise to the difference spectrum) . This is essentially 
a balanced homodyne detection scheme, except that in 
this case the signal arises in the low noise difference chan- 
nel rather than the sum channel. 

At the focus of a high NA lens it is easier to control the 
circular polarisation of a narrowband beam rather than 
a broadband one, and, commonly, the pump and probe 
fields are supplied by the same beam, so E pr = E p . To 
address this we can generalise (j4]) to the case of arbitrary 
ellipticity of both Stokes and pump beams [HJ 

AS cx C(9,cb)I^I s (l - 3p)x JVR Im{xf 111 } (5) 

where 

C{9, <f>) = ^ [sin(40) + sin(20) + sin(20) cos(40)] . (6) 
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Figure 2: (colour online) Interferometric correction of the 
NRB in toluene. S+ and S- are the raw CARS spectra mea- 
sured along the n+ and n_ directions of Fig. [R;b). The 
difference spectrum AS = S+ — S- contains non-dispersive 
lineshapes at the correct Raman shifts (c.f. Fig. [3]). Inset: 
close-up of the region around 1000 cm -1 showing the AS peak 
shifted with respect to the dispersive CARS peaks (lines are 
a guide to the eye). 



9 and cf> characterise the ellipticity of the Stokes and 
pump/probe fields, respectively, and are defined as the 
angle between the fast axis (set parallel to the x-axis) of 
a quarter waveplate and an input linear polarisation. For 
the experimentally practical situation where the broad- 
band Stokes is constrained to be linear (ie C(0,4>)) the 
SIP-CARS signal, AS, is maximised for an elliptical 
pump with (f> = t/8, 3 Vs- Significantly, AS retains the 
same spectral form regardless of the ellipticity of the ex- 
citation beams. In general, the NRB will be removed as 
long as the polarisation ellipses are symmetric with re- 
spect to the measurement axes, while the ellipticity de- 
termines the amplitude of AS. This decoupling of NRB 
removal and signal amplitude simplifies the alignment un- 
der tight focussing conditions. Requirements on polari- 
sation purity are somewhat relaxed and the polarisation 
can be set in situ at the focus by iterativcly minimising 
the NRB and maximising the difference signal at reso- 
nance. 

Experiments with an elliptical pump and broadband 
linear Stokes were performed with the apparatus detailed 
in Figure [ljc). Dispersion in the long internal PCF of 
the source limited temporal overlap between pump/probe 
and Stokes, such that the effective Stokes range was 750- 
2300 cnr 1 (though strong resonances were detectable 
to 3000 cm -1 [39] ) . The pump/probe width was 30 
cm. . This system served as a rigorous demonstration 
of the robustness of the technique, as the PCF generated 
Stokes beam has relatively low spectral coherence, and 
all the wavelengths used were outside the design range 
for the input objective lens which significantly affected 
SNR through poor focussing performance. Single-shot 
broadband interferometric NRB removal for this system 
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would not be possible by any other optical technique. 

Correction of Raman lineshapes and removal of the 
NRB in SIP-CARS is shown for toluene in Figure^ The 
two CARS spectra, S+ and S_, measured at ±45° from 
the Stokes polarisation exhibit the asymmetric disper- 
sive lineshapes and spectrally varying NRB (the varia- 
tion reflecting the spectrum of the Stokes beam) which 
is characterstic of CARS measurements. The difference 
spectrum, AS, shows no NRB and Raman peaks are sym- 
metric and occur at the correct vibrational energy. The 
inset shows this in grea ter detail for the ring breathing 
mode at 1004 cm" 1 |36|. 

Direct comparison to spontaneous Raman spectra 
can be made by equating the depolarisation ratio in 
Equation [4] with the spontaneous Raman depolarisation 
ratio[35l l371| p = R±/Rn, where Ru and R± are sponta- 
neous Raman spectra with incident and scattered polari- 
sation mutually parallel and perpendicular, respectively. 
Then Ru cx Imjxun} and R± cx Im{%i22i}, and from 
Equation HI 



AS 



cx R n - 3R ± . 



(7) 



Note that the normalisation of the difference spec- 
trum by Is corrects the peak amplitudes for variation 
in the Stokes spectrum to allow comparison with the 
spontaneous Raman spectrum; it is not necessary for 
removing the NRB or for quantitative measurements. 
For a non-resonant sample, Xim — Xnr and X1221 = 
XAr_R,/3, so © reduces to P x = xnrE%E s /V2, P y = 

—iXNR,EpE s /3V2 (where we have set E pr = E p ). The 
sum spectrum, YiS. is then 



S5 cx P X P* 



P P* 
r v r v 



-Xhx 



NR.- 



(8) 
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Figure 3: (colour online) Comparison of normalised SIP- 
CARS to the spontaneous Raman spectrum formed by Ru — 
3R± (Ru parallel polarised, R± perpendicular polarised). 
Curves offset for clarity. Top panel: cyclohexane. Bottom 
panel: toluene. The SIP-CARS and spontaneous Raman 
spectra agree closely, to within the resolution of the SIP- 
CARS measurement. Inset: Concentration dependence of 
various resonances of toluene (in methanol: 100% to 6.25% 
toluene) and cyclohexane (in toluene: 100% to 5% cyclohex- 
ane) in the presence of stongly varying NRB (see text). "ESnr 
is T,S integrated over the range 1800-1980 cm -1 . Error bars 
are within marker size. Curves have been offset and scaled for 
clarity (toluene: tl=1004 cm -1 , t2=1208 cm -1 ; cyclohexane: 
cl=1267 cm" 1 , c2: 1444 cm" 1 ). 



The spectral form of Is can therefore be determined from 
that of T,S obtained from a non-resonant medium; we 
have employed a glass coverslip for this purpose as the 
glass resonant response is slowly varying and weakens to- 
wards higher wavenumbers. We find complete agreement 
between SIP-CARS and spontaneous Raman in terms of 
both spectral position and relative peak heights (shown 
for cyclohexane and toluene in Figure [3]). As expected 
from ([5]) , the amplitude of each Raman line is scaled by 
the depolarisation ratio, with peaks going negative for 
resonances with p > providing a powerful approach 
for the differentiation of otherwise similar spectra on the 
basis of mode symmetry. Modes for which p ~ 1 /z do not 
appear, such as the CH3 'umbrella' deformation mode of 
toluene at 1379 cm" 1 . 

As the non-resonant response, xnr, amplifies the reso- 
nant CARS signal, quantitative measurements in hetero- 
geneous media require account to be taken of any vari- 
ation of Xnr] for example, as density changes across a 
sample. SIP-CARS is self-calibrating, in that xnr can be 



monitored in the same measurement if a spectrum con- 
tains a non-resonant region (such as the 'quiet' region 
exhibited by biological samples). From ([8|), measuring 
£5 at a reference frequency away from resonances gives a 
quantity 'SSnr cx Xnr- The SIP-CARS difference spec- 
trum can therefore be normalised by \/SSnr to give a 
signal which is linear in concentration and independent 
of the strength of the non-resonant response. 

We find that the NRB signal from toluene is approx- 
imately twice as strong as for cyclohexane and much 
stronger than for methanol. Therefore, to demonstrate 
linear concentration dependence in a situation with vary- 
ing NRB, a range of mixtures of cyclohexane in toluene 
and toluene in methanol were investigated. For cyclo- 
hexane in toluene the NRB therefore increased with de- 
creasing cyclohexane concentration, while for toluene in 
methanol the NRB had the opposite trend. ESWi?, was 
measured by integrating the sum spectrum over the range 
1800-1980 cm" 1 , and used to normalise the SIP-CARS 
measurement, which displayed a linear signal dependence 
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on concentration (inset of Figure [3]) . 

In summary, by exploiting the third-order polarisation 
response, SIP-CARS allows acquisition of CRS spectra 
free of NRB, with complete agreement to spontaneous 
Raman measurements. Spectra are amplified by the non- 
resonant respone and are quantitative, with a linear con- 
centration dependence. The method uses only passive 
polarisation optics, has low stability requirements, and is 
suitable for any laser system capable of generating CARS, 
permitting single-shot interferometric NRB removal with 
broad PCF generated supercontinua. 

This work was supported by the Biotechnology and 
Biological Sciences Research Council (BB/F016344) and 
the European Metrology Research Programme (NEW02- 
REG2). 
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